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Abstract: We present a discussion of the effects induced by bulk viscosity 
on the very early Universe stabihty. The viscosity coefficient is assumed to be 
related to the energy density p via a power-law of the form C, — (^oP'^ (where 
^0, s = const.) and the behavior of the density contrast in analyzed. 

In particular, we study both Einstein and hydrodynamic equations up to first 
and second order in time in the so-called quasi-isotropic collapsing picture near the 
cosmological singularity. As a result, we get a power-law solution existing only in 
correspondence to a restricted domain of Co- The particular case of pure isotropic 
FRW dynamics is then analyzed and we show how the asymptotic approach to 
the initial singularity admits an unstable collapsing picture. 

1 The Quasi-lsotropic model 

In 1963, E.M. Lifshitz and I.M. Khalatnikov IJ proposed the so-called quasi-isotropic solu- 
tion. This model is based on the idea that a Taylor expansion, in time, of the 3-metric can 
be addressed. In this approach, the Friedmann solution becomes a particular case of a larger 
class of solutions, existing only for space filled with matter. In fact, for an ultra-relativistic 
perfect fluid (p — p/3, p being the energy density of the fluid), the spatial metric assumes 
the form ^ o-ap{x'^) t, asymptotically towards the singularity as t —^ 0. 

As a function of time, the 3-mctric is expandable in powers of t. The quasi-isotropic 
solution is formulated in a synchronous system (i.e., goa = 0, goo = — 1), with a line element 
of the form 

ds^ = ~dt^ + -fap{t, X'^)dx°'dx'^ , 7a/3 = i aa/3 + t^ bal3 + ■■■ , (1) 

and a"^ is defined as a"'^ap^ — (5", moreover 6^ — a^'^bjfj. Using the extrinsic curvature 



kaf3, Einstein Eqs. read [2] 

27f^*(V7«^^)+^f = T^f-K-^f , « =5,lnV7, 

where indicates the Ricci tensor and 7 — det[7Q/3] ~ + 16) det[aQ^]. 

Generalized Line-Element - To include dissipative effects into the evolution of the energy 
source we deal with a more complex (no integer powers) 3- metric expansion 

7a;3 = t"aa/3 + tna/3 + ... , = t-^ a'^P -ty~^-b-P + (2) 

where we impose the constraints for the space contraction (i.e., x > 0) and for the consistence 
of the perturbative scheme (i.e., y > x). 

Viscous corrections - In this work, we treat the immediate generalization of the LK 
scheme considering the presence of dissipative processes within the asymptotic fluid dy- 
namics, as expected at temperatures above 0{lQ^^GeV). This extension is described by an 
additional term in the expression of the perfect fluid energy-momentum tensor: 

Tf,v^ \ P (4u^Ui. + 9tj.v) - C,u'^^p{u^u^, -f 5^^) , uf.p = dtln^ , (3) 

where p = p/3 denotes the usual thermostatic pressure and C, is the bulk viscosity coefficient. 
In particular, here we assume this quantity as a function of the Universe energy density p; 
according to literature developments, we express C as a power-law of the form C, = (0 p", 
where Co is a constant and s is a dimensionless parameter (0 ^ s ^ 1/2) [1]. In what follows, 
we fix the value s = 1/2 in order to deal with the maximum effect that bulk viscosity can 
have without dominating the dynamics [5] , in view of the phenomenological issue of thermal 
equilibrium perturbations which characterizes this kind of viscosity. 

Writing now the 00-components of Einstein Eqs., we can expand the energy density as 
follows 

eo , eib , eib \ 



^2 t2-y+x ' V'- t \ 2eo 

where we have assumed the condition Wq ~ 1 (with uq = —1) whose consistence must be 
verified a posteriori and eo, ei — const.. 

The Energy Density Solution - In order to obtain an analytical expression for the energy 
density, we match the fO-Einstein Eq., after substituting expressions (U), with the hydrody- 
namical ones T^. = 0. It is worth noting that, in the non- viscous case (Co = 0), the energy 
density solution is determined without exploiting the hydrodynamical Eq., since p directly 
comes out from the f 0-gravitational Eq. . In our approximation (ua is neglected with respect 
to Mo), these Eqs. can be combined together and solved order by order in the 1/t expansion 
(in the asymptotic limit t ^ 0). Since for the coherence of the solution we impose y > x, hy 
solving leading-order identities we get 

x=l/[l-2^Co] , eo = |x^ (5) 

The parameter Co has the restriction Co ^ Ysn/s in order to satisfy a; > 0: so the exponent 
X of the metric power-law runs from 1 (which corresponds to Co = 0) to 00. We remark 
that this constraint arises from a zeroth-order analysis and defines the existence of a viscous 
Friedmann-like model. Comparing now the two first-order identities we get 

?/ = 2, ei = -^x^ +2x'^ -2x . (6) 
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We now narrow the validity of the parameter x to the values which satisfy the constraint 
X < y, guaranteeing the consistence of the model. Thus, from ([5]), the quasi-isotropic solution 
exists only if 

Co < Co = 2/3V3 , (7) 

i.e., the viscosity is sufficiently small. For values of the viscous parameter Co that over- 
come the critical one Co i the asymptotic quasi-isotropic solution can not be addressed, since 
perturbations would grow more rapidly than zeroth-order terms. 

We are now able to write the final expression of the energy density and the density contrast 
(6). With the constraint 1 ^ .t < 2, we get 

We note that the density contrast evolution is strongly damped by the presence of dissipative 
effects: this behavior implies that 6 approaches the singularity more weakly as i — > when 
the viscosity runs to Co- In correspondence with this threshold value the density contrast 
remains constant in time and hence it must be excluded by the possible Co choices. 

We conclude this section verifying the consistence of our model. By the analysis of aO- 
gravitational Eq. we get, up to the dominant-order of expansion: Uq. ^ t^~^. The assumption 
Wq ~ 1 is, therefore, well verified: in fact, UaU^ ^ ^6-3x ^^^^ neglected in the contraction 
U/j.u'^ = — 1 and our approximation scheme results to be self-consistent. 



2 The Pure Isotropic FRW Model 

In this Section, we briefly investigate the effects that bulk viscosity has on the stability of the 
pure isotropic Friedmann- Robertson- Walker (FRW) Universe [51 [7] ■ We start here from the 
perturbed Einstein Eqs. in a FRW background filled with ultra-relativistic viscous matter, 
whose coefficient C = Co corresponds to the choice 5 = 1/2. The zeroth-order dynamics is 
described by the energy conservation Eq. and the Friedmann one. Assuming the conformal 
time rj {dt — adrf), the solutions are 

p = Ca-(^+^'^) , a^aiT]^/'^, cj = l-xCo, X = V547rG , (9) 

being C an integration constant and ai = (Scj^ttCG/S)^/^'^. Since we consider an expanding 
Universe, we have w > obtaining the constraint ^ Co < Vx- 

The analysis of perturbations moves from the viscous energy-momentum tensor and the 
scalar representation of the perturbed 5-metric /i^ in a synchronous reference system [2]. 
Solving the perturbed Einstein Eq. we obtain the following expression for the the density 
contrast 5 = ^f/p 

5 ^ [Ci7]3-V^ + + C^ri^-"'"- + C^yi^-"'% (10) 

where Ci,2,3,4 are constants. 

As issue of this analysis, we find that two different dynamical regimes appear when viscosity 
is taken into account and the transition from one regime to the other one takes place when the 
parameter Co overcomes the given threshold value q'^^^^^^ = 1/3^, However, in both these 
stages of evolution, the Universe results to be stable as it expands; the effect of increasing 
viscosity is that the density contrast begins to decrease with increasing time when Co is 
over the threshold. It follows that a real new feature arises with respect to the non-viscous 
analysis when the collapsing point of view is addressed: if Co > Co*'^^^^j the density contrast 
explodes asymptotically and the isotropic Universe results unstable approaching the initial 
singularity. 

We conclude comparing the pure isotropic picture to the quasi-isotropic one since also 
the pure Friedmann-singularity scheme is preserved only if we deal with limited values of 
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the viscous parameter. The threshold values of the the two pictures satisfy the condition 
^*{FRW) ^ ^*/3 (see Eq. ([7]) using geometric units): this constraint is physically motivated 
if we consider, as it is, the FRW model as a particular case of the quasi-isotropic solution. 
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